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INTRODUCTION 


Goutp (1849) DESCRIBED THREE SPECIES of Littorina 
(L. lepida, L. plena and L. scutulata) from the north- 
eastern Pacific coast that have subsequently been con- 
sidered as morphological variants of L. scutulata. On the 
basis of reproductive characteristics, however, MURRAY 
(1979) demonstrated that this taxon was a mixture of two 
species but was unable to identify morphological differ- 
ences between the two species. The confounding of the 
two species has remained a problem, due to the variance 
and degree of overlap in univariate morphometric char- 
acteristics. In this paper I have employed the multivariate 
statistical techniques of discriminant analysis to identify 
morphological differences between the reproductive types 
of Murray (op. cit.) and principal component analysis 
to describe size and shape variation between the two spe- 
cies. It was also possible using discriminant analysis to 
classify existing type material to one of the two species 
and to identify synonymies in Goutp’s (1849) material. 

This study would not have been possible without the 
loan of type material by Dr. J. Rosewater, U.S. National 
Museum of Natural History, Washington, D.C. and 
Dr. K. J. Boss, Museum of Comparative Zoology, Har- 
vard University, Cambridge, Massachusetts for which I 
am grateful. 


METHODS 


Statistical analyses were based on whorl number, shell 
length, whorl height, width perpendicular to the colu- 
mellar axis, maximal width, shell depth and apical angle 
of preserved females of known spawning history from 
Murray’s (1979) study. Linear measurements were taken 
with a pair of vernier calipers to the nearest 0.01 mm (see 


Figure 1). Apical angle (in radians) was derived as the arc 
tangent of the shell length: maximal width ratio. All 
measurements were made with the aperture oriented up- 
wards. The discriminant and principal component an- 
alyses were conducted using the Biomedical Computer 
Programs — P series library version BMDP-77 developed 
by the Health Sciences Computing Facility, UCLA, Los 
Angeles, California. The discriminant analysis (BMDP- 
7M, program revised November 1979) employed a step- 
wise variable entry procedure followed by stepwise vari- 
able removal to maximise the generation of alternative 
models for highlighting morphological differences (see 
DRAPER & SMITH, 1966). The classification of individuals 
to species was based on a jackknife estimation procedure 
in order to minimise bias. Principal components were ex- 
tracted from the covariance matrix of log transformed 
measurements of shell length, whorl height, width per- 
pendicular to the columellar axis, maximal width and shell 
depth after the suggestion of Joticozur (1963) using 
BMDP-4M (program revised November 1979). Whorl 
number could not be included because it is measured in 
units that give it a minimal contribution to the principal 
components (JoLICOEUR, op. cit.). 


Morphological Differences Associated 
with Reproductive Types within the 
Littorina scutulata Complex 


The earlier study of Murray (1979) established that two 
species were being confounded as Littorina scutulata. 
However, without knowledge of the age of individuals the 
differences in the means of the morphometrics could not 
be ascribed to species differences. In addition, the overlap 
in the distributions of each character meant that uni- 
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variate traits could not be used to tell the species apart 
reliably. 


SD 


Figure 1 


Morphological measurements made on each shell. 
AA — apical angle CA ~ columellar axis 
MW - maximal width SD -— shell depth SL — shell length 
W - width WH -— whorl height 
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Discriminant analysis is ideally suited to the resolution 
of such problems. It is an extremely robust procedure that 
allows one to combine morphometrics (e.g., lengths, 
weights, angles, etc.) with morphological traits (e. g., num- 
ber of whorls, presence or absence of a trait, etc.) to con- 
struct divisive criteria amongst groups known to be dis- 
tinct based on some independent criterion. The divisive 
criteria are linear functions of morphological characters 
that yield the widest separation of groups. In the present 
application we know the two species differ in their repro- 
ductive characteristics and hypothesize the existence of an 
unknown set of concomitant morphological differences. 
We do not know a priori what this set will be, hence we 
employ a stepwise variable entry procedure that will pro- 
duce a maximum number of alternative discriminant 
models. We subsequently choose the set of variables that 
yields discriminant functions with the highest proportion 
of correct classifications of reproductive types and of type 
material. If the probability of correctly classifying indi- 
viduals is sufficiently better than by chance, we accept our 
hypothesis that there exist morphological differences be- 
tween the reproductive types. We can then employ these 
functions to identify individuals of unknown reproductive 
history. 

In this study the best discriminant model for separating 
the two species is the following pair of four variable func- 
tions: Z,==15.71X, — 12.21%, + 20.58X; + 4.75%, 
— 51.48; Z,= 13.34X%, — 20.39X. + 33.06X;, + 12.11 
X, — 70.57 where: X, = whorl number, X, = shell 
length, Xs==whorl height, and Xs==shell depth. The 
classificatory power of these functions is indicated in 
Table 1. With these functions 95.6% of the specimens 
known to be Littorina scutulata were correctly classified 
and 96.0% of the specimens known to be different from 
L. scutulata were also correctly classified. Type material 
of L. scutulata, with the exception of 3 of 5 idiotypes 
from San Francisco (MCZ 169360), was also correctly 
classified. It seems probable that this sample and that of 
USNM 796179 are mixtures of the two species, since the 
probability of misclassification is small (4%) and the col- 
lecting locality is an area where both species commonly 
co-occur (MURRAY, 1979). 

Type material of Littorina plena and L. lepida was 
consistently classified as different from L. scutulata. Since 
the type specimens of L. plena and L. lepida cannot be 
distinguished from each other, they should be regarded 
as conspecific. As first reviser I propose the reinstate- 
ment of Littorina plena (to be described later in this 
paper) for the other member of the L. scutulata complex. 
Littorina lepida should be regarded as a synonym of 
L. plena. I regard L. plena as a preferable name, since 


Vol. 24; No. 3 


THE VELIGER 


Page 235 


Table 1 


Classification of material examined based on discriminant analysis of Littorina scutulata 
reproductive types I and II of Murray (1979). 


Number of cases classified as 


% correctly 


Group classified Littorina scutulata 
Littorina scutulata 95.6 43 
Not Littorina scutulata 96.0 2 
MCZ 169222 100.0 0 
MCZ 169289 100.0 0 
MCZ 169360 mixture 21 
USNM 5635 100.0 0 
USNM 5637 100.0 1 
USNM 5640 100.0 1 
USNM 612308 100.0 1 
USNM 677095 100.0 0 
USNM 677096 100.0 0 
USNM 796179 mixture 42 


Not Littorina scutulata 


Comments 


2 Murray’s (1979) Type II? 2 
Murray's (1970) Type 1 $ 9 
Littorina lepida syntypes, Puget Sound 
Littorina plena idiotypes, San Francisco 
Littorina scutulata idiotypes, San Francisco 
Littorina plena type, San Francisco 

Littorina scutulata figured type, Puget Sound 
Littorina scutulata holotype, Puget Sound 
Littorina scutulata paratype, Puget Sound 
Littorina lepida paratype, Puget Sound 
Littorina plena paratype, San Francisco 
Littorina scutulata from San Francisco 


A 
00 


Onroocrwhnd 


Ishell lengths = 10.33, 8.18 mm 
shell lengths = 5.95, 8.55, 8.60, 8.95 mm 


the description of GouLp (1849: 84) more closely describes 
the specimens I have observed than does his description of 
L. lepida (Goutp, 1849: 83) thereby overriding page 
precedence of the Law of Priority (article 24 of the 
ICZN). 

The discriminant functions Zı and Zz have a wider 
utility outside the scope of this paper. Since the two spe- 
cies are easily confused, these functions can be used in 
future studies to improve the accuracy of identification 
using four easily measured shell characters. The following 
example demonstrates how the functions Zı and Zz are 
used. We collect a shell and record its number of whorls 
(X1) as 5.0, shell length (Xz) as 8.10mm, whorl height 
(Xs) as 4.59mm, and shell depth (X4) as 4.69mm. We 
then calculate values for Zı and Zz given these data. 
We classify our specimen according to the following 
inequalities: if Zı < Z the specimen is Littorina scutulata, 
if Zı > Zə the specimen is L. plena. In our example we 
calculate values of Zı= 44.91 and Z2—= 939.51, therefore, 
since Zi > Z2 we predict (with greater than 95% certainty) 
that the specimen is L. plena. This technique is extremely 
powerful and has been used by Murray (1980) to dem- 
onstrate morphological differences within a species due to 
tidal height and parasitism. The reader should bear in 
mind that these functions do not take into account mor- 
phological variation that might exist between sexes, tidal 
level, or differences that result from other environmental 
sources. This cautionary note does not denigrate the utility 


of these functions but they should be used with caution. 
They are included because we expect that differences be- 
tween species are greater than differences within a species. 


Morphological Variation in Size and Shape 
in Littorina scutulata and Littorina plena 


Thus far we have assumed different shell morphologies 
between species and used these differences without char- 
acterizing them. We can now proceed to develop multi- 
variate measures of size and shape for Littorina scutulata 
and L. plena, using principal component analysis. The 
advantage of this technique is that it allows one to develop 
the best single descriptions of mutually independent trends 
in size and shape variability. Best is used in the sense that 
each trend maximally explains the observed variation in 
morphology with respect to that trend. For example, each 
description of size variation (= first principal component) 
maximally describes the variability that can be attributed 
to variation in shell morphology associated with differences 
in size among specimens. The remaining information 
about morphology is explained by variation in shape that 
is independent of size and derived from the residual vari- 
ance of the first principal component. This process is con- 


tinued until the information gained at each extraction 
tails off. 
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Variation in size and shape together explain more than 
98% of the total variation in shell morphology in Littorina 
scutulata and in L. plena, so additional principal compo- 
nents were not extracted. Most of the variation in mor- 
phology is associated with size related variation, as can 
be seen from Table 2 (97.3% for L. scutulata and 97.5% 
for L. plena). Variation in morphology attributable to 
shape is small (1.3% for L. scutulata and 1.9% for L. 
plena). Although shape is a small component of morphol- 
ogy, it is shape variation that exhibits the greatest differ- 
ence between species and probably accounts for our ability 
to separate L. scutulata and L. plena morphologically. 
This becomes readily apparent when we calculate the 
angle between each trend in variation. It should be borne 
in mind that each principal component is a vector (whose 
elements are the coefficients in Table 2) extending through 
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Graphically, we can view each trend as a morpho- 
logical continuum with a greater separation of the species 
along the shape continuum than along the size continuum. 
This separation of species along the shape continuum is 
evident in Figure 2. The ellipses have been drawn in to 
delimit the groups and do not represent confidence ellip- 
ses. This figure does serve to dramatize the different mor- 
phologies for Littorina scutulata and L. plena. We can 
now describe the characteristics of each and highlight the 
differences between species. 


Littorina scutulata Gould, 1849 


holotype: USNM 5640, collected Puget Sound, shell length 
= 12.16 mm 
paratype: USNM 612308, collected Puget Sound, shell 


length = 8.79 mm 


Table 2 


Coefficients of size and shape functions for Littorina scutulata and Littorina plena. 


Littorina scutulata 


size variation 


shape variation 


Littorina plena 


size variation shape variation 


(Ist P.C.) (2nd P.C.) (Ist P.C.) (2nd P.C.) 
Log shell length 0.4907 0.2659 0.5012 0.2152 
Log whorl height 0.4083 —0.2926 0.3740 —0.4460 
Log width 0.4399 —0.4361 0.4216 —0.3823 
Log maximal width 0.4267 —0.3640 0.4187 — 0.3317 
Log shell depth 0.4657 0.7202 0.5062 0.7080 
Variance X 104 229.86 414.11 7.98 
% total variance 97.30 97.53 1.88 


a cluster of points that represents the shell geometry of 
each species. Since it is a vector, we can calculate the 
angle between principal components as a measure of the 
difference between species. If we take Uı:= un, u12, u13, 
u1, us as the vector of coefficients associated with size 
variation in L. scutulata and V, == v11, V12, vis, vis, V15 as 
the vector of coefficients associated with size variation in 
L. plena then the angle between U: and Vı is given by: 
cos == un vu +u: viet. . .+uis vis (JOLICOEUR, 1963). 
The more similar the trends the smaller the angle 0. 
Using the coefficients in Table 2 we calculate an angle 
of 3.1° between the trends in size variation in L. scutulata 
and L. plena. Shape variation between the species is more 
widely divergent and is separated by an angle of 9.9°. 
Shape differs more than size between species. 


Littorina scutulata are small obconic brown snails of 
the mid-littoral zone. Adult snails typically have 5 whorls 
and can be expected to range in size from 11.3mm to 
11.9mm (mean shell length= 11.61 mm). Infrequently 
(15.9% of the snails observed) the base of the body whorl 
will be marked by an indistinct ivory colored band which 
may also be visible as a stripe inside the aperture. The 
apical angle is typically 30.4°. Morphological characters 
are summarized in Table 3. The reproductive features of 
this species definitively identify it (Murray, 1979). Males 
possess a penis with a conspicuous sperm groove running 
dorsally to a sub-terminal bulge. Female L. scutulata pro- 
duce characteristic pelagic egg capsules shaped like in- 
verted saucers 0.7 to 1.0mm in diameter and containing 
I to 14 eggs 105 um in diameter. 
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Size variation 


Figure 2 


Scatterplot of size versus shape variation in Littorina scutulata 
and Littorina plena 


Asterisks (*) represent the mean of each group; X represents Lit- 
torina scutulata from Dillon Beach, California; open circles (O) 
L. scutulata from Rockport, California; closed circles (@) L. plena 
from Newport, Oregon; closed squares (ff) Z. plena from Dillon 
Beach, California. Values of size and shape variation were calcu- 
lated from the coefficients in Table 2 


Littorina plena Gould, 1849 


(synonym: Littorina lepida Gould, 1849) 


holotype: USNM 5635, collected San Francisco, shell 
length = 8.52 mm 

paratypes: USNM 677096, collected San Francisco, shell 
lengths = 6.22, 8.20 mm 


The morphological characteristics of Littorina plena 
are summarized in Table 3. Although this species is mor- 
phologically similar to L. scutulata, adults are generally 
smaller (8.4mm < shell length < 8.9mm as opposed to 
11.3mm < shell length < 11.9mm). Littorina plena as 
noted by Goutp (1849) also tend to be more ovate, having 
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Table 3 


Summary of the morphological characteristics of 
Littorina scutulata and Littorina plena. 
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Littorina scutulata Littorina plena 


mean (+ 1 S.E.) 


Whorl number 5.09 (0.08) 5.40 (0.10) 
Shell length (mm) 11.61 (0.28) 8.60 (0.25) 
Whorl height (mm) 7.00 (0.14) 4.87 (0.11) 
Width (mm) 7.81 (0.17) 5.44 (0.13) 
Maximal width (mm) 6.84 (0.52) 6.35 (0.16) 
Shell depth (mm) 6.37 (0.15) 4.53 (0.14) 
Apical angle (radians) 0.53 (0.04) 0.64 (0.004) 
Size (Ist P.C.) 1.92 (0.03) 1.68 (0.03) 
Shape (2nd P.C.) —0.02 (0.03) —0.18 (0.004) 
Frequency of band on the 
base of the body whorl 15.9% 84.3% 

n 45 50 


an apical angle of 36.7° as opposed to 30.4° in L. scutu- 
lata. The ivory band present but indistinct in 15.97% of 
L. scutulata is more prominent in L. plena and occurs on 
84.3% of the shells studied. The coloration of L. plena 
shells is also more variable than L. scutulata and ranges 
from dark brown to brown tessellated with grey. As in 
L. scutulata, the reproductive features described by Mur- 
RAY (1979) definitively characterize this species. Male 
L. plena have a penis in which the sperm groove runs 
laterally to the tip instead of dorsally as in L. scutulata. 
The penis of L. plena also bears a prominent papilla on 
the dorso-lateral surface proximal to the curvature of the 
penis. This papilla is absent in L. scutulata. Female L. 
plena also differ from L. scutulata in the characteristics 
of their spawn. Female L. plena produce pelagic egg cap- 
sules that resemble automobile wheels slightly greater 
than 1mm in diameter and contain from 4 to 41 eggs 
95.7 wm in diameter. 


CONCLUSIONS 


Discriminant analysis verified the hypothesis that in addi- 
tion to the reproductive differences in Littorina scutulata 
species complex reported by Murray (1979) there exist 
associated morphological differences. Using principal 
components it was possible to demonstrate that most of 
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the difference in morphology between species could be 
attributed to a separation along a multivariate shape con- 
tinuum. However, the primary purpose of this study is to 
provide a way of minimizing the confounding of Littorina 
scutulata and L. plena in future studies. The ease with 
which these species can be distinguished is ultimately 
dependent upon the level of accuracy a researcher is con- 
tent with. If a researcher is content with an approximate 
16% chance of error, then L. plena can be distinguished 
by the presence of an ivory band on the base of the body 
whorl and its generally smaller size. This error rate, how- 
ever, can be reduced to approximately 4% by counting the 
number of body whorls, measuring shell length, whorl 
height and shell depth and employing the discriminant 
functions derived in the text. If a 4% error rate is un- 
acceptable, the two species can be distinguished with 
absolute certainty only by recourse to the reproductive 
characteristics described by Murray (1979). In addition 
to the low error rate, the discriminant functions have the 
additional advantages that they do not require snails in 
breeding condition; they work equally well on empty shells 
as on live material; the necessary measurements are easily 
made and usually of considerable interest in a study; and 
the calculations necessary to determine species identity are 
readily done on a hand calculator. 
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